Universality classes for Coulomb frustrated phase separation 
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Abstract 

We identify two "universality" classes in the Coulomb frustrated phase separation phenomenon. They corre- 
spond to two different kind of electronic compressibility anomalies often encountered in strongly correlated 
electronic systems. We discuss differences and similarities of their corresponding phase diagrams in two- 
and three-dimensional systems. 
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1. Introduction 

A large variety of systems with phase separa- 
tion (PS) tendencies subject to long-range forces 
self-organize in domain patterns [11,13, BH- Re- 
cently, advances in local probe techniques have re- 
vealed mixed states in materials like cuprates and 
manganites 0, 0| rekindling the study of this phe- 
nomenon in strongly correlated electronic systems. 
Indeed it has become clear that strong electron cor- 
relations generally produce a tendency towards PS 
0,i,S,[ifl,[Iii which, however, is frustrated by the 
long-range part of the Coulomb interaction (LRC). 
This leads to the formation of inhomogeneities with 
a typical size determined by the competition among 
long range forces and surface energy effects. 

Since domains have often mesoscopic scales of 
several lattice constant, one can perform a general 
analysis of the frustrated phase separation (EPS) 
mechanism which neglects the microscopic details 
of each specific system while capturing its general 
properties. Tendency towards PS is then recognized 
by the presence of anomalies in the electronic con- 
tribution to the free energy of the system . The 
anomalies found in a large variety of strongly cor- 
related electronic models can be classified in two 



kinds corresponding to a short-range negative com 



la,ll6| or a Dirac- 
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pressibility density region 

[mill 

delta-like negative divergence of the compressibility 
due to the crossing of the free energies of two dis- 



tinct phases [17|, |18[ . Both situations can be cap- 



tured by expanding the short range part of the elec- 
tronic free energy density /e of the system around a 
reference density ric as /e = a\n — nc\'^. Here a < 
encodes the tendency towards PS where 7=1 cor- 
responds to a compressibility divergence and 7 = 2 
to a negative compressibility region. Higher order 
terms are essential to analyze FPS from the limit of 
strong frustration down to the limit of zero frustra- 
tion. As a minimal model we take a contribution 
to the free energy of the form (3{n — n.cf''' that for 
/5 > provide a symmetric double- well form of the 
short range free energy. 



When LRC effects can be considered as a weak 
perturbation upon the ordinary PS mechanism, one 
can achieve a universal picture of the FPS [l^. On 
the contrary, in the strongly frustrated regime, the 
two short-range compressibility anomalies give rise 
to two different behaviors. The aim of this work 
is to review the foremost features of their corre- 
sponding phase diagrams for two- [D = 2) and 
three-dimensional {D = 3) systems embedded in 
the three-dimensional long-range Coulomb interac- 
tion. 
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2. Universal behavior: the weEik frustration 
regime 

In the 7 = 2 case, FPS can be analyzed by means 
of a paradigmatic model augmented with the 
long-range Coulomb interaction. Pattern formation 
within this model (or closely related variants) has 
been considered in a variety of systems [2l[ includ- 
ing mixtures of block copolymers charged col- 
loids in polymeric solutions and electronic sys- 
tems The corresponding Hamiltonian reads: 

« = /"dx[0»(x)-l]' + |V*(x)|" + ^^ 

where the classical scalar field (j> represents the local 
charge density with (p its average density. A rigid 
background ensures charge neutrality. The model 
Eq. ([T]) can be reached by measuring the electronic 
free energy density /e in unit of the barrier height 
a^/ (4/3), densities in units such that the double- well 
minima are dA, (f) — ±1, and distances in unit of the 
bare correlation length ^. This leads to a renor- 
malized Coulomb couphng — 2e^^^~^/(eo|a|) 
with e being the electronic charge and Sq a dielec- 
tric constant due to external degrees of freedom. 

Frustrated phase separation for 7 = 1 is more 
easily described by adding an auxiliary field s, 
equivalent to an Hubbard- Stratonovich variable, 
linearly coupled to the charge. It can be taken as a 
soft [22] or a conventional Ising spin (s = ±1) with 
the sign distinguishing the two competing phases. 
In the remainder we will refer to the latter case 
which is more straightforward to analyze. Domain 
walls of the Ising order parameter are sharp by con- 
struction with a surface tension cr oc 2 J where J in- 
dicates the Ising coupling. Then, in the continuum 
limit, one obtains the following Hamiltonian: 



H ^ J: + J dx[0(x) - s(x)]2 + ^ 

,[0(x)-0] [0(x')-0] 



dx / dx 



(2) 



As before, we measure energy densities in units 
of 0^/(4(3), the minima of the double well are at 
= ±1 and lengths are measured in units of 
^ = Aafi/a^ which represents the analogue of the 
bare correlation length in the present model. More 
precisely it indicates the size that inhomogeneities 
should have for the total interface energy to be of 



the same order as the phase separation energy den- 
sity gain a^/(4/3). As for the model Eq. the 
phase diagram is spanned by the dimensionless av- 
erage density (f> and the renormalized Coulomb cou- 

pling Q2 ^ g2^I3-l/(g^^)^ 

In the absence of the Coulomb interaction, both 
models are subject to ordinary PS in the density 
range \(j>\ < 1 as ruled by the Maxwell construction 
(MC). Hence, the mixed state consists of macro- 
scopic domains with local densities (j) — ±.1. For 
Q 7^ 0, PS is undermined as a thermodynamic phe- 
nomenon since the LRC energy cost grows faster 
than the volume in the thermodynamic limit. Thus 
charged domains at a mesoscopic scale appear. 
Their typical size Id is determined by the compe- 
tition between the LRC cost ~ Q^^d ^^ the 
surface energy density ^ I'l^ ■ These terms are op- 
timized whenever the inhomogeneities get a typical 
size ~ 1/Q2. 

Another important length scale is the screening 
length of the Coulomb interaction that controls the 
relaxation of the electronic charge inside the do- 
mains. Both the models Eqs. ([T]), ([2]), have a char- 
acteristic screening length defined by ~^ ^ 
It thus follows that in the weak frustration regime 
Q << 1, one obtains the hierarchy of length scales 
(in units of ^): >> >> 1 that gives ground for 
a unified treatment of the FPS mechanism. 

The strong separation between the typical size of 
the domains and the typical interface width ^ 1 al- 
lows to consider the smooth interface of the model 
Eq. ID) as sharp with a surface tension defined 
as the excess energy of an isolated interface 23 1. 



Moreover, since the effect of LRC represents a small 
perturbation upon ordinary PS, inhomogeneities 
will appear with local densities near = ±1. Then, 
FPS can be analyzed by expanding quadratically 
the double- well potential of Eq. ([1]) around its min- 
ima thus leading to a practical equivalence among 
the two FPS models. 

A good approximation in the weakly frustrated 
regime lies in assuming a uniform density approxi- 
mation (UDA) in which the local charge density is 
assumed constant [H, [2^, | ' 



21, YM- The result- 



ing phase diagram in this approximation is shown 
with the thick lines at the bottom of Fig. [T] for the 
7 = 2 case in 3D systems. Comparison with exact 
results shows that the UDA gives a very accurate 
description. From Fig. [T] one sees that droplet-like 
domains are the stable morphologies on entering in 
the inhomogeneous region. As 4> approaches the 
origin at fixed Q a topological transition to rod-like 
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Figure 1: The 7 = 2 phase diagram in three-dimensional 
systems. The small dots indicates the Gaussian instability 
line Qg. The thin (thick) lines represent first-order transi- 
tions in the strong (weak) coupling approximation. In the 
two limits they overlap with the corresponding numerically 
determined transition lines from the homogeneous phase to 
droplet inhomogeneities (□), from droplets to rods (A) and 
from rods to layers ('0>). Finally the black circle indicates 
the CP. 



structures and subsequently to layered structures 
occurs. A similar behavior is also expected in 2D 
systems. 

Our computations at weak coupling are varia- 
tional so we can not exclude more complicated 
phases including elongation of domains and "fin- 
gering" as in classical systems Naturally the 
ordered phases will be very sensitive to quenched 
disorder. Also in the absence of quenched disorder 
the ground state may be hard to reach on a cooling 
experiment leading to a glass state. [3] 



3. Universality classes: 
tion regime 



the strong frustra- 



By increasing the renormalized Coulomb cou- 
pling in Eqs. H]), ([2]), inhomogeneous states with 
local densities close to the reference density emerge. 
In this case the behavior of the two FPS models 
is radically different leading to two "universality" 
classes. 

First, we consider the 7 = 2 model Eq. ([T]). We 
restrict to 3D systems but similar ideas applies also 
to the 2D case. By computing the static response 
to an external field in momentum space, we get the 
charge susceptibility at k ^ as: 



X(k) 



27rg2 



2-|-6(/) 



The charge susceptibility has a maximum at the 
characteristic finite wavevector fco — [27r(3^] ^^'^ and 



Figure 2: The 7=1 phase diagram in D=2 (full line) and 
D = 3 (dashed line) for the smectic solution. The central 
region corresponds to the mixed state (MS). The inset shows 
an enlargement at strong frustration near the logarithmic 
singularity of the two-dimensional phase diagram. 



diverges approaching the Gaussian instability line 

2 

[dotted line in Fig. [T] Qg — Qc{l ~ ^icj) ) where 
Qc = 1/ \/27r. This indicates an instability to- 
wards a sinusoidal charge density wave (SCDW) 
with direction chosen by spontaneous symmetry 
breaking. At small Q the Gaussian transition line 
predicts inhomogeneities within the spinodal region 
1^1 < 1/VS. This contrast with ordinary PS at 
Q — which implies a mixed state in the global 
density region |(/)| < 1. The situation has been re- 
cently clarified in Ref. 



21[ where it has been shown 



that inclusion of non-Gaussian terms results in a 
first-order phase transition preempting the second- 
order Gaussian instability line except for the criti- 
cal point (CP) (^, Q) = (0,Qc). The mixed region 
smoothly connects with the macroscopically phase 
separated state. 

Away but close to the CP, the transition is weakly 
first-order with more complicated morphologies. 
Approaching the first-order line from above, inho- 
mogeneities are predicted to form a BCC lattice 
with subsequent topological transitions to a planar 
hexagonal lattice of rods and layered structures [see 
Fig. [T] ] . One then finds the same topology as in 
the weak coupling regime, but now inhomogeneities 
have smooth interfaces in between. They continu- 
ously evolve into sharply defined droplets, rods and 
layers (disks and stripes in 2D) as Q ^ with a 
proliferation of higher order harmonics. 

Next, we analyze the 7=1 universality class. 
Fig. [2] shows the phase diagrams for two- and three- 
dimensional systems respectively for striped and 
layered structures which describe a smectic elec- 



tronic liquid phase that possesses orientational or- 
der and breaks the translational symmetry only in 
one direction. As for the previous model other mor- 
phologies will be competitive away from = 0. In- 
deed at weak coupling, as mentioned above, the two 
phases diagrams are identical. Here, however, we 
are interested in the behavior approaching (f) = 
where the considered morphology is the most sta- 
ble one. 

In 3D systems, one finds a maximum frustra- 
tion degree above which only uniform phases are 
allowed. This is antithetical to the behavior of 2D 
systems where the transition lines diverge logarith- 
mically at = 13, 2^ and thus the system always 
breaks into domains no matter how strong the frus- 
trating effects are. 

The difference between 2D and 3D results can 
be traced back to the different role of screening for 
different dimensionalities. In three dimensions, the 
charge density decays exponentially from the do- 
main interfaces on the scale of the Thomas-Fermi 
screening length whereas in 2D systems it decays as 
a power law. Phase separation energy gain stems 
from the region where the electronic density is sig- 
nificantly different from its average value [24 ■ In 
two-dimensions this is fulfilled everywhere in the 
domains, even far from the boundaries. This al- 
lows for domains with any typical size, even ex- 
ponentially larger than the screening length. On 
the contrary, in three-dimensions the systems gain 
PS energy in a region of width 1^ around the inter- 
face. Regions far from the boundary produce an ex- 
ponentially small energy gain and thus the system 
adjust itself to eliminate them. As a consequence a 
maximum size rule is generally valid that says that 
inhomogeneities cannot have all linear dimensions 
much larger than the three-dimensional screening 
length l^^ . This allows for arbitrary large inhomo- 
geneities in 2D systems since one of the dimension 
is already smaller than the l^^. 

4. Conclusions 

In this work we reviewed the main features of 
frustrated phase separation in charged systems con- 
sidering two kind of short-range compressibility 
anomalies often encountered in strongly correlated 
electronic systems. The effect of long-range forces 
can be measured by introducing a dimensionless 
renormalized Coulomb coupling Q which is a mea- 
sure of the amount of frustration. Frustration tends 
to reduce the range of density where a mixed state 



appears hence stabilizing the homogeneous phase 
at densities where ordinary PS would occur. This 
situation is in accord with thermodynamic measure- 
ments [igf of the uniform two-dimensional electron 
gas. 

In the weak frustration regime, the FPS mecha- 
nism is not sensitive to the particular compressibil- 
ity anomaly and a unified treatment can be reached. 
A series a morphological transitions is generally ex- 
pected resembliiig the situation found in other clas- 
sical systems [J, |3| ■ In this limit the phase diagram 
can be safely determined by means of a simple uni- 
form density approximation [24, 25,, 26, ,22, ,28]. 

On the contrary, at strong frustration, two dif- 
ferent universality classes arise. In systems with a 
negative electronic compressibility region (7 — 2) 
a critical value of the frustration exist Qc for both 
D = 2,3. Close to Qc soft inhomogeneities appear. 
They continuously evolve into the sharply defined 
structures of the weak frustration regime. 

For systems with a cusp singularity in the elec- 
tronic compressibility (7 = 1) the system dimen- 
sionality plays a key role. Indeed a maximum frus- 
tration exists only in three-dimensional systems. 

According to the maximum size rule, domains 
cannot have all linear dimensions much larger than 
the screening length. Therefore mesoscopic do- 
mains are generally expected in systems with small 
compressibility as bad metals, systems close to 
metal-insulator transitions and systems with very 
anisotropic electronic properties as indeed found. 
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